Biomarker Changes Associated With Both Dulaglutide and Cardiovascular Events in the REWIND Randomized Controlled Trial: A Nested Case-Control Post Hoc Analysis.

OBJECTIVE: The glucagon-like peptide-1 receptor agonist dulaglutide reduced MACE in the Researching Cardiovascular Events with a Weekly Incretin in Diabetes (REWIND) trial. This article expores the relationship of selected biomarkers to both dulaglutide and major adverse cardiovascular events (MACE). RESEARCH DESIGN AND METHODS: In this post hoc analysis, stored fasting baseline and 2-year plasma samples from 824 REWIND participants with MACE during follow-up and 845 matched non-MACE participants were analyzed for 2-year changes in 19 protein biomarkers. Two-year changes in 135 metabolites were also analyzed in 600 participants with MACE during follow-up and in 601 matched non-MACE participants. Linear and logistic regression models were used to identify proteins that were associated with both dulaglutide treatment and MACE. Similar models were used to identify metabolites that were associated with both dulaglutide treatment and MACE. RESULTS: Compared with placebo, dulaglutide was associated with a greater reduction or lesser 2-year rise from baseline in N-terminal prohormone of brain natriuretic peptide (NT-proBNP), growth differentiation factor 15 (GDF-15), high-sensitivity C-reactive protein, and a greater 2-year rise in C-peptide. Compared with placebo, dulaglutide was also associated with a greater fall from baseline in 2-hydroxybutyric acid and a greater rise in threonine (P < 0.001). Increases from baseline in two of the proteins (but neither metabolite) were associated with MACE, including NT-proBNP (OR 1.267; 95% CI 1.119, 1.435; P < 0.001) and GDF-15 (OR 1.937; 95% CI 1.424, 2.634; P < 0.001). CONCLUSIONS: Dulaglutide was associated with a reduced 2-year rise from baseline of NT-proBNP and GDF-15. Higher rises of these biomarkers were also associated with MACE.

Clinically important alterations in pharmacogene expression in histologically severe nonalcoholic fatty liver disease.

Polypharmacy is common in patients with nonalcoholic fatty liver disease (NAFLD) and previous reports suggest that NAFLD is associated with altered drug disposition. This study aims to determine if patients with NAFLD are at risk for altered drug response by characterizing changes in hepatic mRNA expression of genes mediating drug disposition (pharmacogenes) across the histological NAFLD severity spectrum. We utilize RNA-seq for 93 liver biopsies with histologically staged NAFLD Activity Score (NAS), fibrosis stage, and steatohepatitis (NASH). We identify 37 significant pharmacogene-NAFLD severity associations including CYP2C19 downregulation. We chose to validate CYP2C19 due to its actionability in drug prescribing. Meta-analysis of 16 independent studies demonstrate that CYP2C19 is significantly downregulated to 46% in NASH, to 58% in high NAS, and to 43% in severe fibrosis. Our data demonstrate the downregulation of CYP2C19 in NAFLD which supports developing personalized medicine approaches for drugs sensitive to metabolism by the CYP2C19 enzyme.

Indicators of Kidney Fibrosis in Patients with Type 2 Diabetes and Chronic Kidney Disease Treated with Dulaglutide.

INTRODUCTION: In the AWARD-7 study in patients with type 2 diabetes and moderate-to-severe chronic kidney disease, once-weekly dulaglutide slowed the decline in estimated glomerular filtration rate (eGFR) and decreased the urine albumin/creatinine ratio compared to insulin glargine at the end of 52 weeks of treatment. In this exploratory post hoc analysis, changes in two fibrosis biomarkers, serum PRO-C6 (type VI collagen formation) and urine C3M (type III collagen degradation), were evaluated. METHODS: In the groups treated with dulaglutide 1.5 mg or insulin glargine (N = 330), serum PRO-C6 and urine C3M were measured using competitive enzyme-linked immunosorbent assays. Biomarker changes were assessed by a mixed-effects model for repeated measures. Pearson correlation analyses were conducted to determine associations between changes in kidney fibrosis biomarkers and eGFR measures at 52 weeks. RESULTS: At weeks 26 and 52 of treatment in the overall population, serum PRO-C6 levels were significantly lower in the dulaglutide group versus insulin glargine group with percent change from baseline of (least squares mean ± standard error) -4.6% ± 1.9 and -0.2% ± 2.2 versus 5.7% ± 2.0 and 8.0% ± 2.3 (p < 0.01), respectively, and urine C3M levels were significantly higher in the dulaglutide group versus insulin glargine group with percent change from baseline of 10.9% ± 8.2 and 20.7% ± 8.8 versus -10.0% ± 6.5 and -16.9% ± 6.4 (p < 0.05), respectively. These findings appeared greater in the subgroup with macroalbuminuria. Serum PRO-C6 negatively correlated with eGFR, while urine C3M positively correlated with eGFR. CONCLUSION: Dulaglutide treatment was associated with biomarker changes that indicated lower type VI collagen formation and higher type III collagen degradation compared to treatment with insulin glargine, suggesting a potential drug effect to reduce kidney fibrosis.

Genomic characterization of intrinsic and acquired resistance to cetuximab in colorectal cancer patients.

Anti-EGFR antibodies are effective in therapies for late-stage colorectal cancer (CRC); however, many tumours are unresponsive or develop resistance. We performed genomic analysis of intrinsic and acquired resistance to anti-EGFR therapy in prospectively collected tumour samples from 25 CRC patients receiving cetuximab (an EGFR inhibitor). Of 25 CRC patients, 13 displayed intrinsic resistance to cetuximab; 12 were intrinsically sensitive. We obtained six re-biopsy samples at acquired resistance from the intrinsically sensitive patients. NCOA4-RET and LMNA-NTRK1 fusions and NRG1 and GNAS amplifications were found in intrinsic-resistant patients. In cetuximab-sensitive patients, we found KRAS K117N and A146T mutations in addition to BRAF V600E, AKT1 E17K, PIK3CA E542K, and FGFR1 or ERBB2 amplifications. The comparison between baseline and acquired-resistant tumours revealed an extreme shift in variant allele frequency of somatic variants, suggesting that cetuximab exposure dramatically selected for rare resistant subclones that were initially undetectable. There was also an increase in epithelial-to-mesenchymal transition at acquired resistance, with a reduction in the immune infiltrate. Furthermore, characterization of an acquired-resistant, patient-derived cell line showed that PI3K/mTOR inhibition could rescue cetuximab resistance. Thus, we uncovered novel genomic alterations that elucidate the mechanisms of sensitivity and resistance to anti-EGFR therapy in metastatic CRC patients.

Aurora A Kinase Inhibition Is Synthetic Lethal with Loss of the RB1 Tumor Suppressor Gene.

Loss-of-function mutations in the retinoblastoma gene RB1 are common in several treatment-refractory cancers such as small-cell lung cancer and triple-negative breast cancer. To identify drugs synthetic lethal with RB1 mutation (RB1 mut), we tested 36 cell-cycle inhibitors using a cancer cell panel profiling approach optimized to discern cytotoxic from cytostatic effects. Inhibitors of the Aurora kinases AURKA and AURKB showed the strongest RB1 association in this assay. LY3295668, an AURKA inhibitor with over 1,000-fold selectivity versus AURKB, is distinguished by minimal toxicity to bone marrow cells at concentrations active against RB1 mut cancer cells and leads to durable regression of RB1 mut tumor xenografts at exposures that are well tolerated in rodents. Genetic suppression screens identified enforcers of the spindle-assembly checkpoint (SAC) as essential for LY3295668 cytotoxicity in RB1-deficient cancers and suggest a model in which a primed SAC creates a unique dependency on AURKA for mitotic exit and survival. SIGNIFICANCE: The identification of a synthetic lethal interaction between RB1 and AURKA inhibition, and the discovery of a drug that can be dosed continuously to achieve uninterrupted inhibition of AURKA kinase activity without myelosuppression, suggest a new approach for the treatment of RB1-deficient malignancies, including patients progressing on CDK4/6 inhibitors.See related commentary by Dick and Li, p. 169.This article is highlighted in the In This Issue feature, p. 151.

Genome-wide RNAseq study of the molecular mechanisms underlying microglia activation in response to pathological tau perturbation in the rTg4510 tau transgenic animal model.

BACKGROUND: Activation of microglia, the resident immune cells of the central nervous system, is a prominent pathological hallmark of Alzheimer's disease (AD). However, the gene expression changes underlying microglia activation in response to tau pathology remain elusive. Furthermore, it is not clear how murine gene expression changes relate to human gene expression networks. METHODS: Microglia cells were isolated from rTg4510 tau transgenic mice and gene expression was profiled using RNA sequencing. Four age groups of mice (2-, 4-, 6-, and 8-months) were analyzed to capture longitudinal gene expression changes that correspond to varying levels of pathology, from minimal tau accumulation to massive neuronal loss. Statistical and system biology approaches were used to analyze the genes and pathways that underlie microglia activation. Differentially expressed genes were compared to human brain co-expression networks. RESULTS: Statistical analysis of RNAseq data indicated that more than 4000 genes were differentially expressed in rTg4510 microglia compared to wild type microglia, with the majority of gene expression changes occurring between 2- and 4-months of age. These genes belong to four major clusters based on their temporal expression pattern. Genes involved in innate immunity were continuously up-regulated, whereas genes involved in the glutamatergic synapse were down-regulated. Up-regulated innate inflammatory pathways included NF-κB signaling, cytokine-cytokine receptor interaction, lysosome, oxidative phosphorylation, and phagosome. NF-κB and cytokine signaling were among the earliest pathways activated, likely driven by the RELA, STAT1 and STAT6 transcription factors. The expression of many AD associated genes such as APOE and TREM2 was also altered in rTg4510 microglia cells. Differentially expressed genes in rTg4510 microglia were enriched in human neurodegenerative disease associated pathways, including Alzheimer's, Parkinson's, and Huntington's diseases, and highly overlapped with the microglia and endothelial modules of human brain transcriptional co-expression networks. CONCLUSION: This study revealed temporal transcriptome alterations in microglia cells in response to pathological tau perturbation and provides insight into the molecular changes underlying microglia activation during tau mediated neurodegeneration.

Genomic Aberrations that Activate D-type Cyclins Are Associated with Enhanced Sensitivity to the CDK4 and CDK6 Inhibitor Abemaciclib.

Most cancers preserve functional retinoblastoma (Rb) and may, therefore, respond to inhibition of D-cyclin-dependent Rb kinases, CDK4 and CDK6. To date, CDK4/6 inhibitors have shown promising clinical activity in breast cancer and lymphomas, but it is not clear which additional Rb-positive cancers might benefit from these agents. No systematic survey to compare relative sensitivities across tumor types and define molecular determinants of response has been described. We report a subset of cancers highly sensitive to CDK4/6 inhibition and characterized by various genomic aberrations known to elevate D-cyclin levels and describe a recurrent CCND1 3'UTR mutation associated with increased expression in endometrial cancer. The results suggest multiple additional classes of cancer that may benefit from CDK4/6-inhibiting drugs such as abemaciclib.

Role of TGF-alpha in the progression of diabetic kidney disease.

Transforming growth factor-alpha (TGFA) has been shown to play a role in experimental chronic kidney disease associated with nephron reduction, while its role in diabetic kidney disease (DKD) is unknown. We show here that intrarenal TGFA mRNA expression, as well as urine and serum TGFA, are increased in human DKD. We used a TGFA neutralizing antibody to determine the role of TGFA in two models of renal disease, the remnant surgical reduction model and the uninephrectomized (uniNx) db/db DKD model. In addition, the contribution of TGFA to DKD progression was examined using an adeno-associated virus approach to increase circulating TGFA in experimental DKD. In vivo blockade of TGFA attenuated kidney disease progression in both nondiabetic 129S6 nephron reduction and Type 2 diabetic uniNx db/db models, whereas overexpression of TGFA in uniNx db/db model accelerated renal disease. Therapeutic activity of the TGFA antibody was enhanced with renin angiotensin system inhibition with further improvement in renal parameters. These findings suggest a pathologic contribution of TGFA in DKD and support the possibility that therapeutic administration of neutralizing antibodies could provide a novel treatment for the disease.

The Rat microRNA body atlas; Evaluation of the microRNA content of rat organs through deep sequencing and characterization of pancreas enriched miRNAs as biomarkers of pancreatic toxicity in the rat and dog.

BACKGROUND: MicroRNAs (miRNA) are ~19-25 nucleotide long RNA molecules that fine tune gene expression through the inhibition of translation or degradation of the mRNA through incorporation into the RNA induced silencing complex (RISC). MicroRNAs are stable in the serum and plasma, are detectable in a wide variety of body fluids, are conserved across veterinary species and humans and are expressed in a tissue specific manner. They can be detected at low concentrations in circulation in animals and humans, generating interest in the utilization of miRNAs as serum and/or plasma based biomarkers of tissue injury. MicroRNA tissue profiling in rodents has been published, but sample an insufficient number of organs of toxicologic interest using microarray or qPCR technologies for miRNA detection. Here we impart an improved rat microRNA body atlas consisting of 21 and 23 tissues of toxicologic interest from male and female Sprague Dawley rats respectively, using Illumina miRNA sequencing. Several of the authors created a dog miRNA body atlas and we collaborated to test miRNAs conserved in rat and dog pancreas in caerulein toxicity studies utilizing both species. RESULTS: A rich data set is presented that more robustly defines the tissue specificity and enrichment profiles of previously published and undiscovered rat miRNAs. We generated 1,927 sequences that mapped to mature miRNAs in rat, mouse and human from miRBase and discovered an additional 1,162 rat miRNAs as compared to the current number of rat miRNAs in miRBase version 21. Tissue specific and enriched miRNAs were identified and a subset of these miRNAs were validated by qPCR for tissue specificity or enrichment. As an example of the power of this approach, we have conducted rat and dog pancreas toxicity studies and examined the levels of some tissue specific and enriched miRNAs conserved between rat and dog in the serum of each species. The studies demonstrate that conserved tissue specific/enriched miRs-216a-5p, 375-3p, 148a-3p, 216b-5p and 141-3p are candidate biomarkers of pancreatic injury in the rat and dog. CONCLUSIONS: A microRNA body atlas for rat and dog was useful in identifying new candidate miRNA biomarkers of organ toxicity in 2 toxicologically relevant species.

In Vivo and In Vitro Characterization of Basal Insulin Peglispro: A Novel Insulin Analog.

The aim of this research was to characterize the in vivo and in vitro properties of basal insulin peglispro (BIL), a new basal insulin, wherein insulin lispro was derivatized through the covalent and site-specific attachment of a 20-kDa polyethylene-glycol (PEG; specifically, methoxy-terminated) moiety to lysine B28. Addition of the PEG moiety increased the hydrodynamic size of the insulin lispro molecule. Studies show there is a prolonged duration of action and a reduction in clearance. Given the different physical properties of BIL, it was also important to assess the metabolic and mitogenic activity of the molecule. Streptozotocin (STZ)-treated diabetic rats were used to study the pharmacokinetic and pharmacodynamic characteristics of BIL. Binding affinity and functional characterization of BIL were compared with those of several therapeutic insulins, insulin AspB10, and insulin-like growth factor 1 (IGF-1). BIL exhibited a markedly longer time to maximum concentration after subcutaneous injection, a greater area under the concentration-time curve, and a longer duration of action in the STZ-treated diabetic rat than insulin lispro. BIL exhibited reduced binding affinity and functional potency as compared with insulin lispro and demonstrated greater selectivity for the human insulin receptor (hIR) as compared with the human insulin-like growth factor 1 receptor. Furthermore, BIL showed a more rapid rate of dephosphorylation following maximal hIR stimulation, and reduced mitogenic potential in an IGF-1 receptor-dominant cellular model. PEGylation of insulin lispro with a 20-kDa PEG moiety at lysine B28 alters the absorption, clearance, distribution, and activity profile receptor, but does not alter its selectivity and full agonist receptor properties.

The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance.

The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R(2)0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.

Transcriptome signature of resistance exercise adaptations: mixed muscle and fiber type specific profiles in young and old adults.

This investigation examined the effects of acute resistance exercise (RE), progressive resistance training (PRT), and age on the human skeletal muscle Transcriptome. Two cohorts of young and old adults [study A: 24 yr, 84 yr (n = 28); study B: 25 yr, 78 yr (n = 36)] were studied. Vastus lateralis biopsies were obtained pre- and 4 h post-RE in conjunction with the 1st and 36th (last) training session as part of a 12-wk PRT program in study A, whereas biopsies were obtained in the basal untrained state in study B. Additionally, the muscle fiber type specific (MHC I and MHC IIa) Transcriptome response to RE was examined in a subset of young and old women from study A. Transcriptome profiling was performed using HG U133 Plus 2.0 Arrays. The main findings were 1) there were 661 genes affected by RE during the 1st and 36th training bout that correlated with gains in muscle size and strength with PRT (termed the Transcriptome signature of resistance exercise adaptations); 2) the RE gene response was most pronounced in fast-twitch (MHC IIa) muscle fibers and provided additional insight into the skeletal muscle biology affected by RE; 3) skeletal muscle of young adults is more responsive to RE at the gene level compared with old adults and age also affected basal level skeletal muscle gene expression. These skeletal muscle Transcriptome findings provide further insight into the molecular basis of sarcopenia and the impact of resistance exercise at the mixed muscle and fiber type specific level.

Reverse-genetic analysis of the two biotin-containing subunit genes of the heteromeric acetyl-coenzyme A carboxylase in Arabidopsis indicates a unidirectional functional redundancy.

The heteromeric acetyl-coenzyme A carboxylase catalyzes the first and committed reaction of de novo fatty acid biosynthesis in plastids. This enzyme is composed of four subunits: biotin carboxyl-carrier protein (BCCP), biotin carboxylase, α-carboxyltransferase, and ß-carboxyltransferase. With the exception of BCCP, single-copy genes encode these subunits in Arabidopsis (Arabidopsis thaliana). Reverse-genetic approaches were used to individually investigate the physiological significance of the two paralogous BCCP-coding genes, CAC1A (At5g16390, codes for BCCP1) and CAC1B (At5g15530, codes for BCCP2). Transfer DNA insertional alleles that completely eliminate the accumulation of BCCP2 have no perceptible effect on plant growth, development, and fatty acid accumulation. In contrast, transfer DNA insertional null allele of the CAC1A gene is embryo lethal and deleteriously affects pollen development and germination. During seed development the effect of the cac1a null allele first becomes apparent at 3-d after flowering, when the synchronous development of the endosperm and embryo is disrupted. Characterization of CAC1A antisense plants showed that reducing BCCP1 accumulation to 35% of wild-type levels, decreases fatty acid accumulation and severely affects normal vegetative plant growth. Detailed expression analysis by a suite of approaches including in situ RNA hybridization, promoter:reporter transgene expression, and quantitative western blotting reveal that the expression of CAC1B is limited to a subset of the CAC1A-expressing tissues, and CAC1B expression levels are only about one-fifth of CAC1A expression levels. Therefore, a likely explanation for the observed unidirectional redundancy between these two paralogous genes is that whereas the BCCP1 protein can compensate for the lack of BCCP2, the absence of BCCP1 cannot be tolerated as BCCP2 levels are not sufficient to support heteromeric acetyl-coenzyme A carboxylase activity at a level that is required for normal growth and development.

Developing and applying a gene functional association network for anti-angiogenic kinase inhibitor activity assessment in an angiogenesis co-culture model.

BACKGROUND: Tumor angiogenesis is a highly regulated process involving intercellular communication as well as the interactions of multiple downstream signal transduction pathways. Disrupting one or even a few angiogenesis pathways is often insufficient to achieve sustained therapeutic benefits due to the complexity of angiogenesis. Targeting multiple angiogenic pathways has been increasingly recognized as a viable strategy. However, translation of the polypharmacology of a given compound to its antiangiogenic efficacy remains a major technical challenge. Developing a global functional association network among angiogenesis-related genes is much needed to facilitate holistic understanding of angiogenesis and to aid the development of more effective anti-angiogenesis therapeutics. RESULTS: We constructed a comprehensive gene functional association network or interactome by transcript profiling an in vitro angiogenesis model, in which human umbilical vein endothelial cells (HUVECs) formed capillary structures when co-cultured with normal human dermal fibroblasts (NHDFs). HUVEC competence and NHDF supportiveness of cord formation were found to be highly cell-passage dependent. An enrichment test of Biological Processes (BP) of differentially expressed genes (DEG) revealed that angiogenesis related BP categories significantly changed with cell passages. Built upon 2012 DEGs identified from two microarray studies, the resulting interactome captured 17226 functional gene associations and displayed characteristics of a scale-free network. The interactome includes the involvement of oncogenes and tumor suppressor genes in angiogenesis. We developed a network walking algorithm to extract connectivity information from the interactome and applied it to simulate the level of network perturbation by three multi-targeted anti-angiogenic kinase inhibitors. Simulated network perturbation correlated with observed anti-angiogenesis activity in a cord formation bioassay. CONCLUSION: We established a comprehensive gene functional association network to model in vitro angiogenesis regulation. The present study provided a proof-of-concept pilot of applying network perturbation analysis to drug phenotypic activity assessment.

Sources of variation in baseline gene expression levels from toxicogenomics study control animals across multiple laboratories.

BACKGROUND: The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variance due to individual, environmental, and technical factors. Meta-analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies could yield useful information on baseline fluctuations in gene expression, although control animal data has not been available on a scale and in a form best served for data-mining. RESULTS: A dataset of control animal microarray expression data was assembled by a working group of the Health and Environmental Sciences Institute's Technical Committee on the Application of Genomics in Mechanism Based Risk Assessment in order to provide a public resource for assessments of variability in baseline gene expression. Data from over 500 Affymetrix microarrays from control rat liver and kidney were collected from 16 different institutions. Thirty-five biological and technical factors were obtained for each animal, describing a wide range of study characteristics, and a subset were evaluated in detail for their contribution to total variability using multivariate statistical and graphical techniques. CONCLUSION: The study factors that emerged as key sources of variability included gender, organ section, strain, and fasting state. These and other study factors were identified as key descriptors that should be included in the minimal information about a toxicogenomics study needed for interpretation of results by an independent source. Genes that are the most and least variable, gender-selective, or altered by fasting were also identified and functionally categorized. Better characterization of gene expression variability in control animals will aid in the design of toxicogenomics studies and in the interpretation of their results.

DNA microarray data integration by ortholog gene analysis reveals potential molecular mechanisms of estrogen-dependent growth of human uterine fibroids.

BACKGROUND: Uterine fibroids or leiomyoma are a common benign smooth muscle tumor. The tumor growth is well known to be estrogen-dependent. However, the molecular mechanisms of its estrogen-dependency is not well understood. METHODS: Differentially expressed genes in human uterine fibroids were either retrieved from published papers or from our own statistical analysis of downloaded array data. Probes for the same genes on different Affymetrix chips were mapped based on probe comparison information provided by Affymetrix. Genes identified by two or three array studies were submitted for ortholog analysis. Human and rat ortholog genes were identified by using ortholog gene databases, HomoloGene and TOGA and were confirmed by synteny analysis with MultiContigView tool in the Ensembl genome browser. RESULTS: By integrated analysis of three recently published DNA microarray studies with human tissue, thirty-eight genes were found to be differentially expressed in the same direction in fibroid compared to adjacent uterine myometrium by at least two research groups. Among these genes, twelve with rat orthologs were identified as estrogen-regulated from our array study investigating uterine expression in ovariectomized rats treated with estrogen. Functional and pathway analyses of the twelve genes suggested multiple molecular mechanisms for estrogen-dependent cell survival and tumor growth. Firstly, estrogen increased expression of the anti-apoptotic PCP4 gene and suppressed the expression of growth inhibitory receptors PTGER3 and TGFBR2. Secondly, estrogen may antagonize PPARgamma signaling, thought to inhibit fibroid growth and survival, at two points in the PPAR pathway: 1) through increased ANXA1 gene expression which can inhibit phospholipase A2 activity and in turn decrease arachidonic acid synthesis, and 2) by decreasing L-PGDS expression which would reduce synthesis of PGJ2, an endogenous ligand for PPARgamma. Lastly, estrogen affects retinoic acid (RA) synthesis and mobilization by regulating expression of CRABP2 and ALDH1A1. RA has been shown to play a significant role in the development of uterine fibroids in an animal model. CONCLUSION: Integrated analysis of multiple array datasets revealed twelve human and rat ortholog genes that were differentially expressed in human uterine fibroids and transcriptionally responsive to estrogen in the rat uterus. Functional and pathway analysis of these genes suggest multiple potential molecular mechanisms for the poorly understood estrogen-dependent growth of uterine fibroids. Fully understanding the exact molecular interactions among these gene products requires further study to validate their roles in uterine fibroids. This work provides new avenues of study which could influence the future direction of therapeutic intervention for the disease.

Underlying mechanisms of pharmacology and toxicity of a novel PPAR agonist revealed using rodent and canine hepatocytes.

Marked species-specific responses to agonists of the peroxisome proliferator-activated alpha receptor (PPAR alpha) have been observed in rats and dogs, two species typically used to assess the potential human risk of pharmaceuticals in development. In this study, we used primary cultured rat and dog hepatocytes to investigate the underlying mechanisms of a novel PPAR alpha and -gamma coagonist, LY465608, relative to fenofibrate, a prototypical PPAR alpha agonist. As expected, rat hepatocytes incubated with these two agonists demonstrated an increase in peroxisome number as evaluated by electron microscopy, whereas the peroxisome number remained unchanged in dog hepatocytes. Biochemical analysis showed that rat hepatocytes responded to PPAR agonists with an induction of both peroxisomal and mitochondrial beta-oxidation (PBox and MBox) activities. Dog hepatocytes treated with both PPAR agonists, however, did not show increased PBox activity but did demonstrate increased MBox activity. Analysis of peroxisomal beta-oxidation gene expression markers by quantitative real-time PCR confirmed that PPAR agonists induced the peroxisomal enzymes, acyl-coenzyme A (CoA) oxidase (Acox), enoyl-CoA hydratase/L-3-hydroxyacyl-CoA dehydrogenase (Ehhadh), and 3-ketoacyl-CoA thiolase (Acaa1) at the transcriptional level in rat hepatocytes, but not dog hepatocytes. Expression of mRNA for the mitochondrial beta-oxidation gene hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase (Hadhb), however, increased in both rat and dog hepatocytes, consistent with biochemical measurements of peroxisomal and mitochondrial beta-oxidation. Repeat-dose nonclinical safety studies of LY465608 revealed abnormities in mitochondrial morphology and evidence of single-cell necrosis following 30 days of dosing exclusively in dogs, but not in rats. Microarray analysis indicated that dog hepatocytes, but not rat hepatocytes, treated with LY465608 had an expression profile consistent with abnormalities in the regulation of cell renewal and death, oxidative stress, and mitochondrial bioenergetics, which may explain the canine-specific toxicity observed in vivo with this compound. This increased sensitivity to mitochondrial toxicity of canine hepatocytes relative to rat hepatocytes identified using gene expression was confirmed using the fluorescent indicator tetramethylrhodamine ethyl ester (TMRE) and flow cytometry. At doses of 0.1 microM LY465608, canine hepatocytes showed a greater shift in fluorescence indicative of mitochondrial damage than observed with rat hepatocytes treated at 10 microM. In summary, using rat and dog primary hepatocytes, we replicated the pharmacologic and toxicologic effects of LY465608 observed in vivo during preclinical development and propose an underlying mechanism for these species-specific effects.

Identification of blood biomarkers of rheumatoid arthritis by transcript profiling of peripheral blood mononuclear cells from the rat collagen-induced arthritis model.

Rheumatoid arthritis (RA) is a chronic debilitating autoimmune disease that results in joint destruction and subsequent loss of function. To better understand its pathogenesis and to facilitate the search for novel RA therapeutics, we profiled the rat model of collagen-induced arthritis (CIA) to discover and characterize blood biomarkers for RA. Peripheral blood mononuclear cells (PBMCs) were purified using a Ficoll gradient at various time points after type II collagen immunization for RNA preparation. Total RNA was processed for a microarray analysis using Affymetrix GeneChip technology. Statistical comparison analyses identified differentially expressed genes that distinguished CIA from control rats. Clustering analyses indicated that gene expression patterns correlated with laboratory indices of disease progression. A set of 28 probe sets showed significant differences in expression between blood from arthritic rats and that from controls at the earliest time after induction, and the difference persisted for the entire time course. Gene Ontology comparison of the present study with previous published murine microarray studies showed conserved Biological Processes during disease induction between the local joint and PBMC responses. Genes known to be involved in autoimmune response and arthritis, such as those encoding Galectin-3, Versican, and Socs3, were identified and validated by quantitative TaqMan RT-PCR analysis using independent blood samples. Finally, immunoblot analysis confirmed that Galectin-3 was secreted over time in plasma as well as in supernatant of cultured tissue synoviocytes of the arthritic rats, which is consistent with disease progression. Our data indicate that gene expression in PBMCs from the CIA model can be utilized to identify candidate blood biomarkers for RA.

Optimization and validation of small quantity RNA profiling for identifying TNF responses in cultured human vascular endothelial cells.

INTRODUCTION: Affymetrix oligonucleotide microarrays are widely used in basic and applied research (Lander, E.S., (1999). Array of hope. Nature Genetics 21, 3-4; Lockhart, D.J. & Winzeler, E.A. (2000) Genomics, gene expression and DNA arrays. Nature 405, 827-836.) The need for a significant amount of starting RNA has limited its use in applications where the amount of RNA is limiting, such as with Laser Captured Microdissection (LCM), small biopsies, or peripheral blood in rodent models. To overcome this limitation, various RNA amplification and labeling methods have been described, however, further optimization and validation of these methods are needed. METHODS: Here we reported using the Arcturus technology to optimize amplification and labeling of small amounts of RNA for Affymetrix microarray studies. We assessed the technical feasibility and variation introduced by differences in starting RNA quantity and differences in technical performance by microarray hybridization. RESULTS: We demonstrated that the current approach is reliable to amplify as little as 40 ng total RNA, and it is suitable for Affymetrix studies yielding satisfactory quantitative chip performance. We also showed that differences in labeling methods contribute more to variation than the differences in starting RNA quantity per se. As a model, we studied the well-documented TNF-induced inflammatory responses in cultured human vascular endothelial cells. We were able to recapitulate the TNF-induced responses using small RNA sample profiling.